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Quantum evaporation of Callan-Giddings-Harvey-Strominger (CGHS) black holes is analyzed in the mean 
field approximation. This semi-classical theory incorporates back reaction. Detailed analytical and numerical 
calculations show that, while some of the assumptions underlying the standard evaporation paradigm are borne 
out, several are not. Furthermore, if the black hole is initially macroscopic, the evaporation process exhibits 
remarkable universal properties (which are distinct from the features observed in the simplified, exactly solu- 
ble models). Although the literature on CGHS black holes is quite rich, these features had escaped previous 
analyses, in part because of lack of required numerical precision, and in part because certain properties and 
symmetries of the model were not fully recognized. Finally, our results provide support for the full quantum 
gravity scenario recently developed by Ashtekar, Taveras and Varadarajan. 
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/. Introduction. Since the early nineties, a number of 2- 
dimensional (2D) black hole models have been studied to gain 
further insight into the quantum dynamics of black hole evap- 
oration. Physically, the most interesting among them is due to 
Callan-Giddings-Harvey-Strominger(CGHS) [1]. Simplified 
versions of this model are exactly soluble but also have im- 
portant limitations discussed, e.g., in [2, 3], Therefore results 
obtained in those models are not reliable indicators of what 
happens in the full CGHS dynamics. In this letter we present 
key results from a new analysis of CGHS black holes using 
a mean-field or semi-classical approximation. These findings 
are surprising in two respects. First, several features of the 
standard CGHS paradigm [2] of quantum evaporation are not 
realized. Second, black holes resulting from a prompt collapse 
of a large Arnowitt-Deser-Misner (ADM) mass exhibit rather 
remarkable behavior: after an initial transient phase, dynam- 
ics of various physically interesting quantities at right future 
null infinity l£ flow to universal curves, independent of the 
details of the initial collapsing matter distribution. This uni- 
versality strongly suggests that information in the collapsing 
matter on 2^ can not in general be recovered at . How- 
ever, we also find strong evidence supporting the scenario of 
[4] in which the 5-matrix from (left past infinity) to 
is unitary. This distinction between unitarity and information 
recovery is a peculiarity of 2D. 

In this letter we summarize the main results. An exten- 
sive treatment can be found in [5]; details about the numerics 
in [6]; and a more thorough investigation of the full quantum 
issues in [7]. 

//. Model. In the CGHS model, geometry is encoded in a 
physical metric g and a dilaton field </>, and coupled to N 
massless scalar fields fa. Since we are in 2D with M 2 topology, 
we can fix a fiducial flat metric rj and write g as g ab = ilr] ab . 
Then it is convenient to describe geometry through $ := e~ 2 * 
and := The model has 2 constants, k with dimen- 

sions [L]^ 1 and G with dimensions [ML] -1 . 

Our investigation is carried out within the mean field ap- 
proximation (MFA) of [4, 7] in which one ignores quantum 
fluctuations of geometry but not of matter. To ensure a suf- 
ficiently large domain of validity, we must have large N and 
we assume that each scalar field /j has the same profile on 1^ . 



singularity..-' '-.last ray 




Figure 1 . A Penrose diagram of an evaporating CGHS black hole in 
the mean field approximation (MFA). The incoming state is the vac- 
uum on XjJ , and left moving matter distribution on . The collapse 
creates a generalized dynamical horizon (GDH), which subsequently 
evaporates. Quantum radiation fills the spacetime to the causal future 
of matter. Inside the GDH, a singularity forms in the geometry. It 
meets the GDH when the latter shrinks to zero area. The "last ray" 
emanating from this meeting point is a future Cauchy horizon. 

Black hole formation and evaporation is described entirely in 
terms of non-linear partial differential equations. Denote by 
z the advanced and retarded null coordinates of r\ so that 
r) a b = 2<9( a z + db)Z~. We will set d± = d/dz . Then we 
have the evolution equations 

□(„)/< =0 «*• □(,)/< = 0. (1) 

for matter fields, and 

d+ cL $ + k 2 Q = G (T+_) = NGhd+ <9_ ln^e- 1 ) 

$<9+ <9_lne = -G(T+_) = -NGhd+d- ln($e" 1 )(2) 

for geometric fields O, <f>. The terms on the right side are 
quantum corrections to the classical equations due to confor- 
mal anomaly and encode the back reaction of quantum radi- 
ation. As in 4D general relativity there are constraints which 
are preserved by the evolution equations: 

-d 2 _ $ + <9_$<9_ ln6 = G(T__ } 

-d% $ + d+ $<9+ In 9 = G (T++) . (3) 
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Here, TV := TV/24 and (T a i,) denotes the expectation value of 
the stress-energy tensor of the TV fields /j. 

We solve this system of equations as follows. As is stan- 
dard in the CGHS literature, we assume that prior to z + = 
the space-time is given by the classical vacuum solution and 
matter falls in from Zp after that (see Fig. 1). Therefore, to 
specify consistent initial data, it suffices to choose a matter 
profile f + (z + ) on Zp , and solve for the initial (0,$) us- 
ing (3). We then evolve (0, $) to the future of the initial data 
surfaces using (2). Trivially, fi(z + , z~) = f + (z + ) from (1). 

We now discuss the interpretation of solutions to these 
equations via horizons, singularities and the Bondi mass. Note 
first that in analogous 4-dimensional (4D) spherically sym- 
metric reductions, $ is related to the radius r by $ = n 2 r 2 
[2, 5]. Therefore, a point in the CGHS space-time (M,g) 
is said to be future marginally trapped if 9+$ vanishes and 
is negative there [2, 8]. The quantum corrected "area" 
of a trapped point is given by a := (<£> — 2NGH). The 
world-line of these marginally trapped points forms a gener- 
alized dynamical horizon (GDH). As time evolves, this area 
shrinks because of quantum radiation (hence 'generalized': 
the world-line is time-like rather than space-like). The area 
finally shrinks to zero. The MFA equations are formally sin- 
gular where $ = 2NGH; thus at the end-point of evaporation 
the GDH meets a space-like singularity. The 'last ray' — the 
null geodesic from this point to — is the future Cauchy 
horizon of the semi-classical space-time. See Fig. 1. 

We assume (and this is borne out by the simulations) that 
the semi-classical space-time is asymptotically fiat at Zp in 
the sense that, as z + — > oo, the field $ has the following 
behavior along z~ = const lines 



$ = A(z~) e KZ + B(z~) + Q{e 



(4) 



where A and B are smooth functions of z~ . A similar ex- 
pansion holds for 9. The physical semi-classical metric g a i, 
admits an asymptotic time translation t a . Its affine parame- 
ter y" is given by e~ Ky = A(z~). Up to an additive con- 
stant, y~ serves as the unique physical time parameter at Zp . 
The MFA equations imply that there is a balance law at Zjt 
[4, 7], motivating new definitions of a Bondi mass M^J^ 
and a manifestly positive energy flux F ATV : 
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2 L dz- 2 Mz- j J 



(5) 
(6) 



so that d{M^J^/dy~ = -F ATV . In the classical theory 
(H = 0), there is no energy flux at Zp, and Mg^^ reduces 
to the standard Bondi mass formula, which includes only the 
first two terms in (5). Previous literature [1, 2, 8-10] on the 
CGHS model used this classical expression also in the semi- 
classical theory. But we will see that this traditionally used 
Bondi mass, M^^j, is physically unsatisfactory. 

///. Scaling and the Planck Regime. It turns out that the mean 
field theory admits a scaling symmetry. To express it explic- 
itly, let us regard all fields as functions of z^. Then, given any 



solution (6, 3?, TV, /+) to all the field equations and a positive 
number A, (A0, A$, ATV, /+) is also a solution (once z~ is 
shifted to z~ + (lnA)//c) [5, 12]. Under this transformation, 
we have 
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(M, F v , a GDH ) -> A(M, F v , a GDH 



where aQDH denotes the area of the GDH, and M is either the 
Bondi mass Mg^j or the ADM mass Madm- This symme- 
try implies that, as far as space-tune geometry and energetics 
are concerned, only the ratios M/N matter. Thus, whether a 
black hole is 'macroscopic' or 'Planck size' depends on the 
ratios M/N and acDH/-^V rather than on the values of M or 
a GDH themselves. Hence we are led to define 

(M*,M£ ondi ,F*) - (M^M^F^/N, 

and m* =M^ ondi | lastray . (7) 

To compare these quantities to the Planck scale we first note 
that there are some subtleties because Gh is dimensionless in 
2D and careful considerations lead us to set Mpj = Hk, 2 /G 
and Tp[ = Gh/ k 2 [5]. We can now regard a black hole as 
macroscopic if its evaporation time is much larger than the 
Planck time. Using the fact that, in the external field approx- 
imation, the energy flux is given by Fjjaw = (NHk 2 /2), this 
condition leads us to say that a black hole is macroscopic if 
M* > Gh Mpi. Note that the relevant quantity is M* rather 
than M, The precise nature of this scaling property was not 
appreciated until recently. For example, in [13] it was noted 
that N could be "scaled out" of the problem and that the re- 
sults are "qualitatively independent of N" whereas in fact for a 
given M they can vary significantly as TV changes. Similarly, 
the condition that a macroscopic black hole should have large 
M/N appears in [9]. But it was arrived at by physical consid- 
erations involving static solutions rather than an exact scaling 
property of full equations. 

IV. Results. Here we describe some key results from numer- 
ical solution of the CGHS equations (l)-(2). We consider 
two families of initial data, most conveniently described in 
a "Kruskal-like" coordinate kx + = e KZ . The first is a col- 
lapsing shell used throughout the CGHS literature so far, 

9 M* 

(df + /dx + y = — s( x +-i/k), (8) 

parameterized by M*. The other is a smooth (f + (x + ) is C 4 ), 
two parameter (TV/*, w) profile defined by 



= 4C ( 1 - e( 
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/to : 



)(x+ - 1/k), 
(9) 

where 9 is the unit step function, w characterizes the width of 
the matter distribution, and M* is related to the ADM mass 
via M* « M*(l + 1.39 w). Unraveling of the unforeseen be- 
havior required high precision numerics [6], which is crucial 
in the macroscopic mass limit that is of primary importance. 
Numerical solutions from both classes of initial data were ob- 
tained for a range of masses M* from 2~ 10 to 16, a range of 
widths from w — to w = 4, and TV varying from 0.5 to 
1000. Since we are primarily interested in black holes which 
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Figure 2. The final mass m* versus the initial mass M* (7) for a 
variety of initial data (8-9). The curve fit to the data is m* = a (1 — 



with a w 0.864, /3 w 1.42, and 7 « 1.15. 



are initially macroscopic, here we will focus on M* > 1 and, 
since the computations did bear out the scaling behavior, on 
the case N = 1. We set h= G=K = 1. 

Our numerical simulations show that, as expected, the 
semi-classical space-time is asymptotically flat at but, in 
contrast to the classical theory, Zp is incomplete, i.e., y~ has 
a finite value at the last ray. However, dynamics also exhibits 
some surprising features which can be summarized as follows. 

First, the traditionally used Bondi mass M^^ s can become 
negative and large even when the GDH is macroscopic. For 
CGHS black holes, negative was known to occur [11] 

but only for black holes which are of Planck size even before 
evaporation begins. For initially macroscopic black holes, the 
standard paradigm assumed that M^j^; is positive and tends 
to zero as the GDH shrinks (so that one can attach a 'flat cor- 
ner of Minkowski space to the future of the last ray). Second, 



while the improved Bondi mass, M 



ATV 
Bondi 



does remain positive 



throughout evolution, at the last ray it can be large. In fact 
this 'end state' exhibits a universality shown in Fig 2 where 
to*, the final value of (Aij*, ondi ), is plotted against the rescaled 
ADM mass M* for a range of initial data. It is clear from the 
plot that there is a qualitative difference between M* > 4 and 
M* < 4. In the first case the value of the end point Bondi 
mass is universal, m* w 0.864. For M* < 4 on the other 
hand, the value of to* depends sensitively on M*. Thus in the 
MFA it is natural to regard CGHS black holes with M* > 4 
as macroscopic, and those with M* < 4 as microscopic. Nu- 
merical studies have been used in the past to clarify proper- 
ties of the CGHS model [3, 10, 11, 13], such as dynamics of 
the GDH. However, they could not uncover universal behavior 
because, in the present terminology, they covered only micro- 
scopic cases (M* < 2.5 in all prior studies). 

Third, for macroscopic (A/* > 4) black holes that form 
promptly, after some early transient behavior, dynamics of 
physical quantities at the GDH and at Zj£ approach univer- 
sal curves. By promptly, we mean the characteristic width of 
the ingoing pulse is less than that of the initial GDH (more 
precisely, w/M* < 0.1). This is most clearly demonstrated 
in the behavior of the flux F*, or equivalently the Bondi mass 
Afg ondi , measured at X^. An appropriately shifted affine pa- 
rameter «7 = y~ + const provides an invariantly defined 
time coordinate and Fig. 3 shows the universality of evolution 
of F* and A/g ondi with respect to it. The shift aligns the y~ 




Figure 3. F* and M*; ondi of Eq (7) plotted against j/~T for solutions 
with several values of parameters M* and w of Eqs (8)-(9). In all 
cases F* starts at in the distant past (fty^ <C —1), and then joins 
a universal curve at a time that depends on the initial mass. The 
time when the dynamical horizon first forms is marked on each flux 
curve (which is later for larger w, though note that the mass and 
flux curves for all the M* = 6 cases are indistinguishable in the 
figure). We have not yet found an extrapolation of the flux to the 
last ray — that conclusively answers whether it is finite there. 
However, all functions we tried that fit the data well have a finite 
integrated flux. Moreover, when the flux starts rising rapidly, we are 
still well within the regime where the numerical solution converges, 
and we can follow the solution clearly into a regime where the mass 
has reached its final value of m* w 0.864. 



coordinates amongst the solutions, which we are free to do 
as y~ is only uniquely defined to within a (physically irrele- 
vant) additive constant. Finally, we note that this universality 
is qualitatively different from the known uniqueness results 
for solutions of certain simplified soluble models [14]. It oc- 
curs only if the black hole is initially macroscopic, formed by 
a prompt collapse. And in this case, after the transient phase, 
the behavior of physical quantities at Zp does not depend even 
on the mass. 

The overall situation with universality bares some paral- 
lels to the discovery of critical phenomena at the threshold 
of gravitational collapse in classical general relativity [15] 
where universal properties were discovered in a system that, 
at the time, seemed to have been already explored exhaus- 
tively. Of course, numerical investigations cannot prove uni- 
versality; here we only studied two families of initial data. 
However, since these families, in particular the distribution, 
are not 'special' in any way, we believe this is strong evidence 
that universality is a feature of the 'pure' quantum decay of a 
GDH, pure in the sense that the decay is not contaminated by 
a continued infall from Z^ . 

Finally, along the last ray, our simulations showed that cur- 
vature remains finite. Thus, contrary to wide spread belief, 
based in part on [3], and in contrast to simplified and soluble 
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models, there is no 'thunderbolt singularity' in the metric. 

V. Conclusions. In the external field approximation, the en- 
ergy flux is initially zero and, after the transient phase, quickly 
asymptotes to the Hawking value Fn a w = Nhn 2 /2 = 0.5 
for the constants used in the simulations shown here. In the 
MFA calculation, on the other hand, at the end of the tran- 
sient phase the energy flux is higher than this value, keeps 
monotonically increasing and is about 70% greater than Fn aw 
when Msondi ~ 2NMp\ (see Fig 3). One might first think 
that the increase is because, as in 4D, the black hole gets hot- 
ter as it evaporates. This is not so: For CGHS black holes, 
Taaw = nh/2-K and k is an absolute constant. Rather, the 
departure from Fnaw = 0.5 shows that, once the back reac- 
tion is included, the flux fails to be thermal at the late stage of 
evaporation, even while the black hole is macroscopic . This 
removes a widely quoted obstacle against the possibility that 
the outgoing quantum state is pure in the full theory. 

In the classical solution, is complete and its causal past 
covers only a part of space-time; there is an event horizon. 
But Zp is smaller than X^ in a precise sense: z~, the affine 
parameter along is finite at the future end of X^. This is 
why pure states on l£" of a test quantum field /_ on the clas- 
sical solution evolve to mixed states onX^ [4, 7], i.e., why the 
S matrix is non-unitary. In the MFA, by contrast, our analysis 
shows that as expected y~ is finite at the last ray on X^. Thus, 
1^ is incomplete whence we cannot even ask if the semi- 
classical space-time admits an event horizon; what forms and 
evaporates is, rather, the GDH. However, this incompleteness 
also opens the possibility that 2"^, the right null infinity of 
the full quantum space-time, may be larger than Zj£ and uni- 
tarity may be restored. Indeed, since there is no thunderbolt, 
space-time can be continued beyond the last ray. In the mean 
field theory the extension is ambiguous. But it is reasonable to 
expect that the ambiguities will be removed by full quantum 
gravity [16]. Indeed, since we only have (0.864/24)Afpi of 
Bondi mass left over at the last ray per evaporation channel 
(i.e., per scalar field), it is reasonable to assume that this re- 
mainder will quickly evaporate after the last ray and MgJ^i 



and F ATV will continue to be zero along the quantum exten- 
sion Zj£ of 1^ . Form of F ATV now implies that X^ is ' as long 
as' as IjJ and hence the S'-matrix is unitary: The vacuum state 
on evolves to a many-particle state with finite norm on X^ 
[4, 7]. Thus unitarity of the S matrix follows from rather mild 
assumptions on what transpires beyond the last ray. 

Note, however, this unitarity of the S'-matrix from X^ to the 
extended X^ does not imply that all the information in the in- 
falling matter on X^ is imprinted in the outgoing state on T£ . 
Indeed, the outgoing quantum state is completely determined 
by the function y~ (z~) and our universality results imply that, 
on , this function only depends on Madm an d not on fur- 
ther details of the matter profile [5]. Since only a tiny fraction 
of Planck mass is radiated per channel in the portion of 2^ 
that is not already in X^, it seems highly unlikely that the re- 
maining information can be encoded in the functional form 
of y~(z~) in that portion. Thus, information in the matter 
profile on X^ will not all be recovered at even in the full 
quantum theory of the CGHS model. This contradicts a gen- 
eral belief; indeed, because the importance of y~(z~) was not 
appreciated and its universality was not even suspected, there 
have been attempts at constructing mechanisms for recovery 
of this information [9]. 

To summarize, in 2D there are two distinct issues: i) uni- 
tarity of the S -matrix from X£ to ; and ii) recovery of the 
infalling information on X^ at Zj. The distinction arises be- 
cause right and left pieces of I ± do not talk to each other. 
In 4D, by contrast, we only have one X" and only one X + . 
Therefore if the S-matrix from X~ to X + is unitary, all infor- 
mation in the ingoing state at X~ is automatically recovered in 
the outgoing state at X + . To the extent that the CGHS analysis 
provides guidance for the 4D case, it suggests that unitarity of 
the S-matrix should continue to hold also in 4D [7]. 
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